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ABSTRACT 

The data from the last seven experiments performed on polarized deep-inelastic 
scattering on proton and neutron (or deuteron) targets have been analyzed in search of 
a precise determination of the spin fraction carried by the quarks in the nucleon. We 
find that this fraction can be of the size expected from naive quark model arguments, 
provided the gluon axial anomaly is explicitly included and the isosinglet axial charge 
normalization is fixed at a suitably low momentum scale, such that a) the running, 
strong coupling constant is about unit, and b) the orbital angular momentum inside 
the nucleon vanishes. 

We also find that, despite the appeal of this solution of the "nucleon spin 
crisis" , a solution where the axial anomaly is absent and its effects are traded for an 
appreciable strange quark polarization can not however be excluded — because of 
the limited accuracy of the data — unless this latter and/or the gluon polarization 
in the nucleon are explicitly measured. 

1. Introduction. 

Just after the 1988 publication by the European Muon Collaboration (EMC) 
at CERN of their preliminary^ results on the asymmetry in muon polarized deep- 
inelastic scattering (PDIS) on a polarized hydrogen (actually spin-frozen ammonia) 
target, the particle physics community was somewhat shocked to learn that the total 
sum of the quark spins in the proton seemed to be very close to zero, rather than 
somewhere between 1/2 and 3/4, as generally expected from naive quark-model 
arguments. 

As for any unexpected discovery in our field, to begin perhaps with that of 
the muon, the following years have seen both the planning and running of new 
experiments, and the deepening of the theoretical studies (not completely free of 
much heated controversies) trying to clarify this mystery, often know under the 
name of "nucleon spin crisis" . 

The aim of this paper, relying heavily on the results thus accumulated, is 
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to use i) all PDIS data taken on both proton and neutron (or deuteron) targets, 
ii) current phenomenological ideas on the behaviour of the parton distributions at 
both large and small values of the Bjorken variable x, and iii) the results from 
perturbative QCD (PQCD) at next-to-next-to-leading order (NNLO), to produce 
an internally consistent estimate of the quark spin content of the nucleon. 

It will turn out that one can not, at the present level of experimental accuracy, 
unambiguously separate the different flavour components from PDIS data alone, 
nor directly verify the PQCD predictions for the isoscalar sum rule: namely, one 
can not decide on a purely experimental basis on the nature of its unitary-singlet 
component, given the (basically two) different possibilities offered in principle by 
the stage at which one decides to send the quark masses (used as regulators in the 
calculations of the splitting functions) to zero. 

Despite this persistent lack of experimental proof of all PQCD arguments, it 
can however be shown that, for the "most natural" (in the parton model framework) 
of these two possibilities, one can interpret the data as being consistent with a 
negligible intrinsic strange component of the quarks' spins at a suitably low mass 
scale, such that as — 1, consistent with the naive quark model expectation. 

There are two primary reasons why this result contradicts the initial findings 
from the EMC data^: the first, of experimental nature, is that it is now clear 
that the values given by the EMC were too low because of their choice of the 
data employed to normalize them; the second, of theoretical nature, is that 
PQCD corrections were included only at leading order, while it is only at next-to- 
leading order that peculiar features of the polarized deep-inelastic scattering (and 
particularly for its unitary-singlet piece) become to emerge, as it will be clear in 
next section. 

The rest of this paper will be divided in three parts: a presentation of the 
PQCD NNLO corrections to the first-moment isovector and isoscalar sum rules, a 
careful enumeration of the constraints that can be imposed on the various polarized 
distribution functions of the quarks, in terms of which one can describe the polarized 
structure functions Qi'"', and then a short discussion on the results obtained fitting 
to the data simple parametrizations satisfying these constraints. This discussion 
will be centred on the spin composition of the nucleon, and in particular on the 
need (or absence of any need) for an intrinsic strange component As in it. 

2. The First Moments of t/f " and PQCD at NNLO. 

The sum rules for the first moments Iq'^{Q^) = Jq dxgi'^{x,Q^) are often 
discussed separating the non-singlet, {p — n) combination from the isoscalar one, 
{p + n), which contains both a non-singlet and a singlet part. Actually, despite 
the technical difficulties involved in dealing with the latter, both sum rules have 
equal footing in PQCD and should not be considered as separate entities, apart 
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from historical, or practical, considerations. 

Their only difference lies indeed in the fact that the axial coupling involved 
in the first is extremely well measured from neutron /3-decay, while the couplings 
involved in the second are outside direct measurement, and can be arrived at only 
through more or less founded theoretical arguments. It is therefore expedient to 
separate them, since the first can either be considered a good test of higher-order 
PQCD, or, alternatively (and this will be the attitude taken in the present paper), 
a useful normalization for the non-singlet part of the moments /o'"(<5^). 

The first-moment sum rule for the difference between proton and neutron po- 
larized structure functions, known as the Bjorken sum rule (BjSR), reads^ (written 
in its full QCD garb) 

lE^'iQ') = dx[g^,{x, Q') - g^x, Q')] = ^ ■ Csia,) ■ Qa + {h.t.),,^, , (1) 

where the coefficient C%{as) includes the PQCD corrections to the parton-model 
result, and the symbols (/i.t.)/t=o,i stand for the higher-twist contributions^ (HTC), 
behaving as integer inverse powers of . Since experiments extend from < >= 
2.0 GeV^ (SLAC E142 Collaboration^) to the 10.7 GeV^ of the EMC^'^ inclusion 
of these effects turns out to be crucial to a test of BjSR, Eq. (1). This is indeed 
also the attitude taken by Ellis and Karliner^ in their analisys of PDIS data. HTC 
are also essential^ to connect the PDIS sum rules for proton and neutron targets to 
the Gerasimov-Drell-Hearn sum rules^ at Q"^ = 0. 

An estimate of the HTC was originally given — and revised — by Balitskii, 
Braun and Kolesnichenko"^ , and then repeated by several groups of authors^ '^"^r in 
this paper preference will be given to the estimate by Ross and Roberts^, mostly 
because of its widespread use in other analyses of the same data. 

Let us now focus our attention on some aspects of the recently improved "'^"'^ 
evaluation of C^{as) and on their implications. This coefficient is now known to 
better than 0(q;^), and to NNLO accuracy it can be written as 

Cs{as) = 1 - ^ • [1 + cf\N,) . ^ + cf\N,) . + ...] , (2) 

TT TT TT 

where the numerical values of cf'l are listed, for Nf from 3 to 5, in Table I, and their 
complete expressions can be found in the original paper by Larin and Vcrmaseren^^. 
Due to the large values of the constants cf'l for Nf — 3, 4, Cgicts) is evidently 
decreasing with ag much more steeply than expected on the basis of the leading 
order estimate Cg = 1 — ttg/yr at low values of Q^, such as those of SLAC experiments 
E1424 (< g2 2.0 GeV2) and EU3^^'^^ (< 3 GeV^), and of the 

preliminary deuteron data from the Spin Muon Collaboration^^ (SMC) at CERN 
(< (52 4 g GeV2). 
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An initial comment is in order here: actually, the expansion in ag for C8{as) is 
known to one order beyond^^'^^ what has been written in Eq. (2), but its unitary- 
singlet partner, Ci{as), has been estimated^^ only to 0{a^), as most of the quanti- 
ties computed in PQCD to extract the coupling ctgi for its consistency, any PQCD 
analysis must be performed to the same, fixed order in the coupling, using the 
/^-function at that order to express the running of the coupling as{Q'^). 

Accordingly, here only NNLO expansions will be used, together with the ex- 
pansion for a.siQ'^) (and the estimate of its scale Aj^{Nf)) at the same order: a 
recent PQCD analysis by Bethke and Catani^^ of all available data (both space- 
and time-like, but with no inclusion of PDIS ones) led to a conservative estimate 
of A^(5) = 200 ± 50 MeV, and therefore to the value A^(3) = 410 ± 100 MeV 
that will be employed throughout. There is a reason to overlook more recent com- 
pilations: since these tend to include also ctg determinations from PDIS data to 
enrich the available points at low momenta, one should exclude these latter deter- 
minations from the analysis^^ to be free of possible internal biases, and to treat 
PQCD corrections as external inputs rather than as parameters to be fit to the 
data. This is simply and easily accomplished, rather than by re-doing all fits, by 
using the Ajjg{3) from the slightly older compilation of Ref. 17. 



Table I 

Coefficients of higher QCD corrections to BjSR 





ciiNf) 


C2iNf) 


3 


3.5833 


20.2153 


4 


3.2500 


13.8503 


5 


2.9167 


7.8402 



A question better addressed at this point is the actual value of Nf, the number 
of active flavours, to be used at each in the PQCD expansions of the coefficients 
Ci,8{cts) and of the anomalous dimension of the unitary-siglet axial charge. For 
time-like Q^, there is no ambiguity, since for heavy quarks the flavour thresholds 
at = 4m| can be fixed by setting the mass of the i-th. flavour quark ap- 
proximately equal to that of its lowest-mass pseudoscalar meson. In deep-inelastic 
scattering, however, a flavour is active only when appreciably contributing to the 
moment sum rules, i.e. when produced a) in a really inclusive manner (read: not 
only in low-multiplicity events), and b) over an appreciable range in Bjorken's vari- 
able X (say up to a; ~ 1/3). If one sets the beginning of the scaling region at ~ 
2 GeV^ (as indicated by the "classic" SLAC-MIT experiments), the previous re- 
quirements ask for a ~ 17 GeV^ for charm to be an active flavour in DIS, and 
so mandates Nf = 3 for the PQCD analysis of all available PDIS data. 

The choice of Nf = 3 (rather than 4) has no great effect on Cgicts), as one 
can read from Table I, but for the Q^-evolution of the unitary-singlet piece the 
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coefficient of the first-order term in the expansion in powers of as of its integrated 
anomalous dimension almost cancels, for Nf = 4, the first-order one in Ci{as). It 
is therefore important, for the extraction of the quarks' spin fraction, 

Nf 

S = ^ Agj = Au + Ad + As for iV/ = 3 , (3) 

i 

to use the value of Nf appropriate to the range of values of Q"^ where the actual 
data have been taken. This, unfortunately, has not always been done consistently 
by some authors^^. 

In 1974, Ellis and Jaffe^*^, faced with the problem of how to use a sum rule 
akin to Eq. (1) with only hydrogen data, used parton-model ideas, fiavour SU(3) 
symmetry and the Okubo-Zweig (OZI) rule to derive a sum rule for the first moment 
of g{ alone; as stated above, the PQCD-corrected version of such a sum rule (when 
freed of OZI~rule restrictions) is a part of PQCD as fundamental as the BjSR. For 
the isoscalar combination of PDIS structure functions this sum rule becomes the 
(PQCD corrected) EUis-Jaffe sum rule (EJSR) 



CiQ') = f dx[gl{x, Q') + g^x, Q^)] 
Jo 



= ^ ■ Csias) -gs + l' Ci(a,) • goiQ^) + {h.t.)i^=o , (4) 

not to be confused with the original EJSR^°, which was derived for Iq{Q'^) only, 
and without any PQCD corrections to the parton-model plus OZI-rule predictions. 

Additional complications with respect to the BjSR, Eq. (1), arise from the 
facts a) that the isoscalar axial charges of the nucleon are not directly measurable, 
and Qs is indeed derived via flavour-symmetry arguments (apart from symmetry- 
breaking effects), and b) that the unitary-singlet one go{Q'^) couples not only to 
the quarks, but also to the gluons via the axial anomaly and possesses therefore 
anomalous dimensions^ so that its evolution with is not exausted by the PQCD 
coefficient Ci{as) and must be explicitly computed. 

Since this coupling is scheme dependent, this point has been the focus of a 
very heated theoretical debate^^. To cut a long history short, one can summarise 
it by saying that, in conventional parton language where the masses of the partons 
rrii are neglected with respect to the momentum scale Q, and wishing to identify S 
with the spin fraction carried by the quarks actually present in the target proton, 
one has to put 

go{Q') = J:-Nf^AG{Q') , (5) 

where AG is the first moment of the gluon polarized distribution function SG{x) = 
G+{x) — G-.{x), and determine its evolution via the equation (where t — logQ^///^) 

■^9o{t) = -Nf^'ygg{as)go{t) , (6) 
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which relates to the anomalous dimension of the axial anomaly via 



-Nf^^gqM = IggM - /5(Q!s)^ (6') 

where the various •jij = {g,q}) represent the coefficients giving dE/dt and 

dAG/dt in terms of E and AG, whose matrix is diagonalized (in any scheme where 
quark mass regulators are sent to zero) building the combination in Eq. (5) on one 
side, which evolves anomalously as in Eq. (6), and leaving S on the other, free of 
anomalous dimensions since the constraint in Eq. (6') makes the determinant of 
the 2x2 matrix of the (final) coefficients (whose eigenvalues are the anomalous 
dimensions of the two operators mixed by the evolution) vanish identically. This 
has been verified step by step at next-to-leading order^^: that it should hold also 
at higher orders is inferred from the fact that conserved (and partially conserved) 
charges should be free of anomalous dimensions on one side, and on the other 
there is nothing, apart from the mixing with the unitary singlet of the pure gluonic 
world, to distinguish an SU(3)- from a U(3)-symmetric fermionic world, so that 
it is "natural" to expect all purely fermionic operators to be free of anomalous 
dimensions. 

After integrating in as from a normalization scale fi^ to Q^, one obtains the 
integrated anomalous dimension at NNLO 

5'o(/^ ) 33-2iV/ TT 

■^^^24^ 36 8(153 - 19iV/)^ tt 

with the NNLO calculation by Larin^^ of the anomalous dimension. 

As one can read from the above expression, fi'o(Q^) can be drastically reduced 
(still at Nf = 3) from its value at the normalization scale for ^ fi^; on the 
other hand, one can not set //^ — oo and thus drop a.sifi'^) from Eq. (7) for two 
reasons: first, that the anomaly contribution to Eq. (5) is not definable in this 
limit^^'^^, due to AG{Q^) diverging asymptotically as as{Q^)~^ from Eqs. (5) and 
(7), and, second, that we can not keep Nf = 3 while letting — > oo. 

There is another point which deserves consideration: the expression in Eq. (7) 
is of second order in ctg, while it was said from the beginning that the analy- 
sis of the PDIS data will be done consistently at NNLO, or 0{a^). Indeed the 
anomalous dimension (a^) has been computed to 0{ag) by Larin^^: it is how- 
ever a logarithmic derivative in the variable t, which to be integrated has to be 
divided by the /3-function to make the variable change from t to ctg, losing in 
this step a power in the coupling as- This gives no problem if the anomalous di- 
mension is left where it belongs, i.e. as the argument of an exponential function. 
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9oiQ^) = ^o(a*^) ■ 6xp[7(q!s(Q^)) — 7(q:s(/x^))], where 7(0;^) is the integrated anoma- 
lous dimension, but has caused unnecessary problems whenever exp7(Q;s) has been 
expanded in powers of the coupling, since then one had either to go one step further 
in the expansion of the anomalous dimension to find the "missing" power^^'^^, or 
to truncate "prematurely" the series in the coeflBcient Ci(as). Both procedures are 
inconsistent from a PQCD point of view, since they tend to mix one order with the 
next, while the running coupling ag must be defined order by order: it is clear that 
in this way one will end up using ctg at a given order in at least one perturbative 
expansion calculated at a different order (not to be confused with the power of 
as appearing in the Gnal expression, which could depend, as is the case here, on 
additional mathematical manipulations). 

Another, equally unnecessary, but luckily only semantic problem has been 
created by the persons who first misnamed the scale /U^ in the same equation a 
renormalization scale, while it is clear enough that one has to do with just a normal- 
ization scale, which has absolutely nothing to share with the actual renormalization 
procedure^^. The scale //^ is thus free both to appear explicitly in the physical 
expressions, and to be chosen to follow the author's (or authors') theoretical preju- 
dices; certain precautions must however be followed: for instance, working at a fixed 
number of flavours Nf = 3, it would be rather unwise to set it to inflnity, where Nf 
will be six at the best of our knowledge, so that the go{oo)^f^^ so obtained would 
hardly be connected to the physical limit of go{Q^) for — * 00. Note also that 
this go{oo)^f^^ (whatever its meaning) can not be identified with E as defined by 
Eq. (3): if one believed the singlet axial coupling go{Q'^) in the EJSR, Eq. (4), 
to be given by E alone, from the diagonalization of the unitary-singlet operators 
one had also to put 7 = accordingly. In line of principle, these two definitions 
for goiQ"^) could be distinguished on the basis of the different evolutions with 
they predict for the EJSR integrals Iq~^^{Q'^): making the two to coincide at a scale 
Qq, the difference for the EJSR between the case of a running go{Q^) and that of 
gfQ = S at a scale would correspond to 

AC^iQ') = ^ . C^{a,{Q')) ■ {eMlMQ')) - iMQl))] - 1} • 9o{Ql) ■ (8) 

We shall return later on this point, when analysing our numerical results in Section 
4. 

The evaluation of the singlet coefficient Ci{as) runs along the same lines as 
that of the non-singlet one C8{as) and it can be similarly expanded in powers of 
as- 

TT TT TT 

There are hovewer additional complications, for this coefficient receives contribu- 
tions from graphs, not included in the other, whose number grows with the order of 
the calculation. Thus Kataev^^ has produced only an estimate of the constant c^2^ 
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for Nf = 3 (note that one of the conditions under which this was possible breaks 
down for Nf = 4). 

To close this section, the numerical values of the constants for the perturbative 
expansions in orders of (dg/Tr) of Cs{as), Ci{as), and 7(cts) = Sfe7fe • (cks/t'')'^, are 
listed in Table II. Two points must be noted before going to the next section: the 
constants in the singlet part are systematically lower than in the corresponding, 
non-singlet one, giving a slower evolution with Q^, and the trends of the inte- 
grated anomalous dimension exp7(Q;s) and of the coefficient Ci{as) run in opposite 
directions, tending to some extent to compensate each other. 

Also, the first line in Table II makes clear enough the essential difference be- 
tween leading- and higher-order PQCD treatments: besides introducing an anoma- 
lous dimension for go{Q'^), the latter break strongly the accidental, lowest-order 
degeneracy of the coefficient functions Ci^s{o(s)- 



Table II 

Constants in the perturbative expansions at NNLO 



k 




^k 


Ik 





1 


1 





1 


3.5833 


1.0959 


0.6667 


2 


20.2153 


~ 3.7 


0.9907 



3. Measurements and parametrizations for gi{x). 

What is actually measured by experiments in PDIS are not the structure 
functions gi themselves, but rather the polarization asymmetries Ai, related to the 
polarized cross sections a'^\a'^^ by 

D-A^ = (aTT - a^i)/{a^^ + a^i) , (10) 

where D is the target polarization fraction, and g2 is neglected: gi is then related 
to Ai by 

It is obvious that the factor 2x in the denominator makes a direct determi- 
nation of (7i at X — > impossible for finite-accuracy data. The evaluation of the 
EJSR, Eq. (4), depends thus on the parametrization assumed to extrapolate gi 
to X = 0: the usual treatment of this point has till now assumed it to extrapolate 
smoothly to a constant, as gi a + (3 ■ x, in accord with the pion-pole trajectory 
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intercept being close to zero. However, since one does not expect the sea distri- 
butions to couple dominantly to an isovector, pseudoscalar trajectory, but rather 
to an isoscalar one such as the eta, one should rather have for these a behaviour 
x~°'^^^\ with aniO) ~ —1/4, which, together with the expected negative sign for 
the sea contribution, produces a spike in the isoscalar part of as a; — > 0, of the 
type gi ~ a — ^-x^^^ (with a, /3 > 0), perhaps just appearing^^ at very low values of 
the Bjorken variable x in the proton data taken by the Spin Muon Collaboration at 
CERN at < >= 10.0 GeV^. Of course, such a spiky behaviour does not show in 
the integrand of the BjSR, which can therefore be extrapolated smoothly to x = 
according to the conventional practice in this matter. 

This point could have a non-negligible influence on the evaluation of the EJSR, 
Eq. (4) , raising the low-a; contribution to its left-hand side well above the conven- 
tional estimates. A special comment is in order here on the EMC published data^ 
for g^: rather than using them, one should use instead only their values for with 
an adequate set of values for F2 and Rp. Indeed, the EMC calculated g^ using 
a) the ratio predicted by PQCD, systematically smaller than experiment since 
finite-mass corrections are dominant at low values of Q'^ and x (though the effect of 
this choice is not too important at < >— 10.7 GeV^), and b) tiieir values for the 
unpolarized structure function F2, systematically lower than those by the BCDMS 
collaboration^^ (and than the recent New Muon Collaboration (NMC) data^^ as 
well) by as much as 13% at the lowest values of x. Even using their measured 
values of Ai together with a phenomenological parametrization for (a;,Q^) (and 
RP{x, Q^)) to produce g^ix, Q^) at a reference, fixed value of (and assuming Ai 
to vary little^'^^'^^ with Q^, an assumption which a recent analysis from the E143 
Collaboration at SLAC seem to corroborate^^) is not completely free of the above, 
last source of error: indeed only the latest, post-NMC parametrizations^^ have 
dropped the unpolarized EMC data altogether, while all previous analyses ended 
up averaging over the two, conflicting sets of data for small values of x. We have 
therefore re-normalized the published EMC data^ for gf with the known ratio of 
BCDMS to EMC data, and this will be the data set referred to as EMC^^^ in the 
rest of the paper. 

As just mentioned two paragraphs above, to build the first moments the actual 
measurements of g^'^ have to be extrapolated in a; to a; = and to a; = 1, to cover 
parts of the integration range not coverable by the experiments on the asymmetry 
A^'^. The second extrapolation does not pose any problem, since Ai tends to 
unit (and R to zero) in the limit x ^ 1 (this is a well-known feature"^'^ of the 
nucleon wavefunction: at high values of x the proton — neutron — structure is 
dominated by the u+ — d+ — quark distribution: some parametrizations^-*^ violate 
this constraint gaining thus additional but unphysical freedom in dealing with the 
small sea components), and Eq. (11) is thus reducing simply to gi ~ F2/{2x), 
whose behaviour as x — ^ 1 is largely determined by well known "counting rules" of 
conventional parton models. 
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In this paper, instead of extrapolating separately in the two extreme ranges of 
X, the choice has been made to fit the data to simple functional forms incorporating 
at least three elements: a) the counting rules, b) reasonable Regge (or other, QCD- 
motivated) behaviours for x — > (separately for valence and sea contributions), 
and c) constraints on the integrals of the parton distribution functions coming from 
a connection between the constituent-quark picture and the quark-parton model, 
originally introduced by Altarelli, Cabibbo, Maiani and Petronzio^^, and recently 
recovered in this context by Fritzsch'^"^ . 

In this picture, the structure functions gff are decomposed in terms of the 
helicity distribution fuctions Sq^ for each active quark flavour (qi = u, d, s), where 
the prime indicates that the contribution from the axial anomaly, formally of order 
ttg, has been included in each flavour's sea distribution, dq[ = 6qi + kgg (E) 6G (the 
symbol (S> will stand from here on for a convolution integral in Bjorkcn variables), 
with kqg the appropriate gluon-to-quark splitting function. It has been sometimes 
said in the literature that the anomaly contributes to gi only at very small x values, 
so that it should be hardly seen in the a;-regions covered by experiments: this is 
true (in the scheme where quark masses rrii are sent to zero) only if the polarized 
gluon distribution is assumed to peak at a; = 0, as e.g. in the intrinsic gluon 
distribution proposed by Brodsky and Schmidt^^. Unfortunately, their distribution 
for 5G{x) has the wrong Regge behaviour for a; ~ 0, contrary to their statement, 
since it requires dominance of SG by the pion trajectory, with intercept 0177(0) ~ 0, 
while it should be dominated instead by a pseudoscalar-glueball trajectory, with 
an expected intercept adO) ~ 077/(0) — —1. When this constraint is imposed on 
the Brodsky-Schmidt formulae, ceteris paribus, most of the anomaly contribution 
falls in the x interval covered by the experiments (e.g. 80% of it in the case of 
the EMC a;-range), making it virtually indistinguishable from the other, intrinsic 
sea distributions, while the integral AG remains of the same magnitude as in Ref. 
34, being solidly tied to the momentum fraction < xq >— 1/2 carried by the 
gluons at the momentum scale < >slac— 1 ~ 2 GeV^ at which these intrinsic 
components are deflned. 

Separating further valence (uy, dy) from sea {u'g, d'g, s') components one can 
write, with self-explanatory notations. 



+4 • [5dy{x,Q'') + 5d',{x,Q'') + 5s'{x,Q'')] , 



(12) 



and 



+4 • [^uy{x, Q2) + Su',{x, Q2) + 5s' {x, g2)] . 



(13) 
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The valence and sea distributions (with the latter primed to distinguish them from 
the intrinsic ones, free of the gluonic contribution, which relate to the quarks spin 
contents /S.qi entering e.g. in the angular momentum sum rule), will in general be 
expressed by the general functional form Sq'y^si.^^Q'^) — ' (1 ~ x)'^^'" ■ 

Pj{x, Q^), where the sum runs over the Regge singularities assumed to dominate the 
distribution at x ~ 0, and the Pj{x, Q^) will be polynomials in x with Q^-dependent 
coefRcients^^'^^. 

For the powers n„ s the parton-model counting rules^^ give n = 2Ns — 1 
(where Ng is the number of "spectator" quarks), or n^; = 3 and Ug — 7 (for gluons 
the rule gives Ug = 5: their further suppression in the distributions Sq'g at high x- 
values comes from the splitting function): it is hovewer known^^'^^ that the valence 
u-quark dominates at high x values over the d-quark in their unpolarized distribu- 
tions, and this fact is commonly "explained" as a consequence of the Pauli principle, 
barring two quarks of the same quantum numbers from being close to each other 
in phase space, and often expressed as a "penalty factor" (1 — a;) in the odd flavour 
distribution function. The same mechanism should operate in the polarized valence 
distributions as well, as in the sea distributions (both polarized and unpolarized), 
suppressing here u-quarks with respect to d-quarks, the "penalty" being now paid 
when an u-antiquark is produced at x ~ 1; while there is undisputable evidence of 
this effect for unpolarized valence distributions in the ratio , tending almost 

linearly to the value 1/4 as a; —> 1, the same effect in the unpolarized sea ones could 
be responsible^^ for the defect of the so-called Gottfried sum rule (one can simply 
check that the figures are indeed of the right order of magnitude, though a detailed 
model would require a complete re-fitting of all unpolarized parton distributions). 
For the PDIS data, this Pauli-principle "penalty factor" will be included only in 
the polarized vaience distributions, for the effects of its presence in the sea ones 
would be so small vis-a-vis the experimental errors that its inclusion would only 
lead to unnecessary mathematical complications in the fitting procedures. 

With this additional factor omitted, and reducing the polynomials Pj{x, Q'^) to 
Q^-dependent factors, independent of x, the sea contributions to the PDIS structure 
functions will reduce to an isoscalar term, simply expressed as 

gn^, Q')sea = P{Q') -x'^-il- xY , (14) 

where we have assumed all sea components to couple dominantly to the r^-meson 
trajectory with intecept a,j(0) ~ —1/4. Note that under this hypothesis the sea 
would contribute to the BjSR only through the difference between u-quark and 
d-quark seas eventually brought about by the Pauli-principle "penalty factor": 
once this latter is dropped, the BjSR contains only the valence distributions and is 
therefore offering a mean to define their normalizations in terms of g^- 

This is accomplished by using the connection between constituent-quark and 
parton model languages of Refs. 32-33, i.e. by writing du = dU ® Suu + SD ® 6ud 
(and so on for all active fiavours and for gluons as well), where the notation Sqq rep- 
resent the g-fiavoured parton polarized density inside the Q-fiavoured constituent 
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quark. Assuming SU(2) symmetry to hold for the partons inside the constituent 
quarks [Q — U, D only), and integrating over all Bjorken variables, one can put 
Auu = Ado = A{qQ)y + A{qQ)s, Adu = Aud = A{qQ>)s and Asu = Asd = As, 
and one finds the relation for the spin content of the valence partons 

^ = ^=-4 (15) 
Ady AD ^ ' 

from the constituent quark model results AU = 4/3, AD = —1/3: note that, since 
these constituent quarks are structured objects, these values do not imply gA = 5/3 
[modulo small recoil corrections), as in naive treatments of the constituent quark 
model. The same hypothesis leads only to the generous bounds —1/4 < Au'^/ Ad'^ < 
1 on the non-strange sea distributions, from the expected sea spin-component ratio 
< A((Jq)s/A((Jq')s < 1 for sea partons inside a constituent quark, where the upper 
and lower bounds are justified respectively by the Pauli-principle "penalty factor" 
mentioned above and by the similarity in their production mechanisms. 

With the constraint of Eq. (15) imposed on the polarized valence distributions 
joined with the neglet of the small isovector part in the polarized sea ones, the 
distributions Su^ and Sd^ can be normalized to the BjSR, independent of their 
functional forms, since Eq. (1) can be reduced to 



/^-"(q2) = dx{6uy - 5dy) = ^ ■ Csias) ■ gA + {h.t. 



For these forms, three parametrizations will be adopted, to check the systematic 
effects on the EJSR integrals of the behaviour assumed in the isoscalar valence 
distributions as x ^ 0. 

The first parametrization (which will be labeled FRP, for fully Reggeized 
parametrization) breaks down each of the polarized distributions for the valence 
quarks into an isoscalar and isovector part, the first dominated at x ~ by the 
?7-meson trajectory, and the second by the pion one, so that, after introducing the 
Pauli-principle "penalty factor" (1 — x) in Sdy and imposing the conditions in Eqs. 
(13) and (14), with the ^-polynomials Pj{x,Q^) again reduced to Q^-dependent 
coefficients, one gets 

(x, Q') = CiQ') . [|1 . 5(^,4)-^ . + |] . (1 - ,;)3 , (17) 

Sdi'Hx, Q^) = cm . [|1 ■ B{\, A)-' . - . (1 - x)4 , (170 

where the normalization C{Q'^) = Cs{ois) ■ Qa + Q' {h.t.)i^i is fixed from Eq. (16) to 
automatically satisfy the BjSR, and the function B{a,P) is the well known Euler's 
beta function. 
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For the second set of model distribution functions (to be labeled SRP, or 
simplified Regge-pole parametrization) one takes the limit aT,(0) — > Q!7r(0) — > 0, 
still keeping the same constraints, so that Eqs. (17), (17') become 

1 fi 

8u^^\x,Q^)^^-C{Q^)-{l-xf , (18) 
o 

5d^:f\x,Q^) = -C{Q^)-{l-x)\ (18') 

where the normalization is again the same. This parametrization has the advantage 
over the previous one of producing an isoscalar structure function + non- 
vanishing in the limit a; — > 0, but at the cost of violating the isotopic spin properties 
of the Reggeon exchanges. 

A third set of model distribution functions has been built following the sug- 
gestion by Bass and Landshoff"^^ that a non-perturbative, two-gluon exchange term 
might develop at small x a behaviour gi ~ log(l/a;), obviously only in the isoscalar 
part, for the PDIS structure function gi. Thus one can write, again with the same 
number of parameters and the same constraints, the functions (labeled, from the 
names of the above-mentioned authors, as BLP) 

54^) (X, Q^) = cm . log ^ + 1^) • (1 - (19) 

54^) (X, Q^) = cm . log i - gil) . (1 - 0.)^ (19') 

this time the r^-meson Regge-pole term has been left out to avoid including one more 

parameter and so keep all three parametrizations on equal statistical footing. Taking 
the PQCD coefficients listed above (with the value for Kj^{h) of Ref. 17) and the 
HTC of Ref. 9 as inputs, there is only one free parameter left for each set of data, 
i.e. the sea normalization P{Q'^). Note that, to account for possible systematics 
either in the data or in the theoretical inputs (particularly the possible inadequacy 
both of the PQCD approximation and of the phcnomenological parametrization 
used, as well as of the evaluation of the HTC), different parameters will be used for 
different data, even when these latter have been normalized at the same value of 
<Q2>. 



4. The isosinglet sum rule and the quark spins in the nucleon. 

Before turning to fitting the parametrizations to the data available (as of 
November 1995) on the PDIS structure functions gi^ it is better to consider the 
information we possess of the axial couplings appearing in the right-hand sides of 
the first-moment sum rules, Eqs. (1) and (4). All pieces of information available 
have been summarized in Table III, which deserves some comments. Usually most 
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analyses^^ of these couplings stop at the information coming from asymmetries in 
the decay products momenta and/or polarizations: while this is enough in most 
of the cases, in some no such measurements are either available or possible, and 
disregard of the information coming from the rates can have drastic effects, even 
on the correctness of the conclusions inferred. Enough to say that dropping the 
information on the AS' = 0, S — > A couplings leads to the conclusion that flavour 
SU(3) symmetry breaking effects are negligible"^^ , relying in fact on a single point, 
being the coupling of the S to both A and E not as well known as those of the latter 
two to the nucleon. Accordingly, Table III includes evidence from both asymmetries 
and rates, while the detailed analysis of these latter can be found elsewhere*^. One 
further thing evident from Table III is that the two sources of information yield fully 
compatible values for the axial couplings, contrary to the statement of Jaffe and 
Manohar'^^, which originated from a completely outdated treatment of the decay 
rate data. 

Table III 

Data on octet baryon axial couplings 

Decays from rates^'^ from asymmetries^^ SU(3) formulae^^ 

n pe-Ve 1.2553 ± 0.0018 1.2573 ± 0.0028 F + D 

S-^Ae-I7e 0.723 ± 0.017 D - ^/StancpF 

E+ Ae+Ue 0.750 ± 0.094 D + VStan^F 

E- ^ ne-Ve - 0.330 ± 0.023 - 0.340 ± 0.017 F-D 

E" ^ nfj,-V^ - 0.237 ± 0.078 F-D 

K^pe-V^ 0.729 ± 0.011 0.718 ± 0.015 F + i(l + ^ tan(/))L» 

A^p^-V^ 0.756 ± 0.139 F + i(l + ^ tan 

H-^Ae-I7e 0.265 ± 0.044 0.250 ± 0.050 F - ^(1 + ^ tan(^)L» 

E-^AiJ-V^ 0.76 F- 1(1 + ^ tan 

E- ^E°e-Ve 1.216 ± 0.147 F + (1 - ^ tan 



Table III does not contain, for reasons of internal clarity, an additional piece 
of information, coming from the study of the asymmetries in the E~ — > A /3-decay, 
i.e. the ratio^^ dv/dA = — 0.01 ± 0.10, which can be considered a bound on the 
size of the E'^ — A mixing effects, which have been included following the analysis 
by Karl'^'^, and fixing tan(/) to the value given by him. 

The SU(3)-symmetric fit to all data listed in Table III, plus the ratio Qv/qa 
for the decay E~ Ae~z7g, yields the parameters F = 0.4678 and D = 0.7876, or 
^8 = 3F - D = 0.616 ± 0.022 a.nd gA = F + D = 1.2554 ± 0.0020 (note that this 
latter acts almost as a constraint due to the very high precision of neutron data 
on both asymmetries and rates); the of the fit is not very high with respect to 
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the number of data points, being 20.17 (about 2 units better than without S°-A 
mixing) versus 15, but is concentrated almost exclusively in the E~ — * A rate, 
whose datum on the axial coupling lies 3a below the SU(3)-symmetric fit value. 
Coming this datum from a good-quality experiment^^, this fact has to be taken 
as initial evidence for some flavour SU(3) breaking in these couplings, being the 
mixing required to explain this discrepancy more than five times the one calculated 
by Karl^^, and thus giving a ratio Qv/qa more than twice its experimental la limit. 
This must be remembered when using fiavour symmetry to extract the nucleon 
spin composition. A further point raised^^ in this context is also the possible effect 
(particularly on the AS' = 1 transitions) of the tensor components in the weak 
axial current: the good agreement between the two columns of data in Table III 
seems to indicate that the effect is not as large as indicated by some fits^^, since 
it would affect the axial couplings extracted from the asimmetries differently from 
those extracted from the rates. 

Since it is the aim of the present paper to concentrate on the analysis of 
the spin composition of the nucleon, the PDIS structure functions will not be fit 
treating both C{Q'^) and P{Q^) as free parameters, but rather taking the first 
from the right-hand side of the BjSR, Eq. (1), and fitting only the second to the 
data, experiment by experiment. Table IV will present the expectations for the 
right-hand side of the BjSR in NNLO PQCD, including the HTC of Ross and 
Roberts^, together with the integrals evaluated by two experimental groups, the 
E143 Collaboration at SLAC^^ and the Spin Muon Collaboration at CERN"^*^. In 
the same table, we shall also display the two expectations for the right-hand side 
of the EJSR obtained assuming validity of the OZI rule, i.e. S = gfg, and either a) 
the naive parton model identification go = T, (i.e. decoupling the contribution from 
the axial anomaly, as is the case for massive quarks), or b) the definition in Eq. (5) 
with the normalizations scale //^ fixed so that a.sip'^) = 1 (and < >= 0): here 
gs is assumed equal to the above fiavour SU(3) symmetry prediction. The choice 
a) amounts to using an isoscalar, PQCD corrected version of the original EJSR^o, 
while b) serves to give an idea of the effect induced by a very reasonably sized axial 
anomaly contribution^^'^^. 

From this table one can read two facts, already mentioned above: the running 
with of the EJSR is slower than that of the BjSR, and therefore harder to see 
in the data unless higher precisions are reached, and the difference between the 
"anomalous" and "non-anomalous" versions of the EJSR, described by Eq. (8) 
(apart from the lower asymptotic value of the first, which can always be traded for 
a breaking of the OZI rule, i.e. a As ^ 0, in the second), is even smaller, and thus 
much harder to see. 
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Table IV 

Sum rule expectations for BjSR and naive EJSR 





EJSR (a) 


EJSR (b) 


BjSR 


BjSR integrals 


Ref. 


2.0 


0.1278 


0.0917 


0.1492 






3.0 


0.1378 


0.0993 


0.1636 


0.163 ± 0.019 


13 


4.6 


0.1448 


0.1045 


0.1733 






6.0 


0.1479 


0.1068 


0.1776 






8.0 


0.1506 


0.1088 


0.1813 






10.0 


0.1523 


0.1099 


0.1836 


0.194 ± 0.038 


46 


10.7 


0.1528 


0.1103 


0.1843 






12.0 


0.1535 


0.1108 


0.1853 







In Table V there will be displayed the parameters P{Q^) obtained from the 
fits to the seven sets of PDIS data analyzed here: the re-normalized EMC set of 
data on proton^'^ at < Q'^ >= 10.7 GeV^, the E142 neutron data"^ at < Q"^ >= 
2.0 GeV2, the SMC prehminary deuteron data^^ at < >^ 4.5 GeV^, the SMC 
proton data^"^ at < Q'^ >— 10.0 GeV^, the E143 data on proton^^ and deuteron^^ 
at < Q2 >= 3.0 GeV^, and the SMC new deuteron data^° at < >= 10.0 GeV^. 
The very preliminary data of the SLAG- Yale E80 and E130 Collaborations^^ have 
not been included in the fits, since they cover either very low values of or a very 
limited range in x values, and would have been, besides, of very little statistical 
significance; also, the recently appeared E143 data^^ are not included, since they 
have become available only during the final redaction of the present paper. The 
factor 1/2 • (1 — fo'D), with ojd = 0.05, has been explicitly included for deuteron 
data, so that all P{Q'^) obtained from the fits have the same normalization. All 
are relatively large and negative, implying a strong reduction from the pure va- 
lence contribution to the EJSR, which under our constraint Eq. (15) would be 
lo^^'iQ^ai = If{Q^)vai = I • C{Q^), i.e. equal to the BjSR integral: from 
this and the previous table one can see that some negative sea is already expected 
even at the level of a naive formulation of the (PQCD corrected) EJSR. The table 
displays also systematic variations in the parameter both with the model used for 
the valence distributions and with the nature of the target, variations which in the 
opinion of the author cast some doubt on the possibility of really testing PQCD in 
the isoscalar combination: since the two kind of variations are of the same order 
of magnitude, and comparable to the statistical errors from the fits, their origin is 
difficult to trace. 
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Table V 

Values of P{Q'^) from fits to the data 

(shown in parenthesis are the x^^values per data point) 



Expt. 


QVGeV^ 


FRP fit 


SRP fit 


BLP fit 




2.0 


-0.38 ±0.07 (0.7) 


-0.35 ±0.07 (0.5) 


-0.23 ±0.13 (3.1) 


E143P 


3.0 


-0.69 ±0.04 (1.9) 


-0.72 ±0.05 (2.2) 


-0.73 ±0.08 (6.2) 


E143'^ 


3.0 


-0.53 ±0.04 (1.3) 


-0.54 ±0.04 (1.2) 


-0.51 ±0.08 (4.7) 


SMC^ 


4.6 


-0.53 ±0.24 (0.6) 


-0.59 ±0.24 (0.7) 


-0.66 ±0.27 (1.1) 


SMC^' 


10.0 


-0.70 ±0.18 (1.7) 


-0.75 ±0.14 (1.0) 


-0.80 ±0.13 (0.7) 


SMC^ 


10.0 


-0.44 ±0.18 (2.7) 


-0.49 ±0.18 (2.5) 


-0.53 ±0.22 (3.8) 


EMC^e. 


10.7 


-0.74 ±0.18 (0.4) 


-0.79 ±0.18 (0.8) 


-0.82 ±0.20 (1.1) 






149.7 / 102 


152.3 / 102 


373.7 / 102 



These results for the sea parameter P{Q^) can now be turned into values for 
the EJSR integrals, since these can be written as 

= \ ■ C{Q^) ± 2 • S(^, 8) • P(g2) . (20) 

Note however that the right-hand side of this equation corresponds to what is 
actually interpolated by the fits only in the case of a deuterium target, when 
^oiQ'^) = I ■ (1 ~ ^'^d) ■ Iq^^{Q'^)' for a proton or neutron target the fits measure, 
respectively, 

= I ■ C{Q') + B{^, 8) • P(Q^) (200 

in the proton case, and 

I^{Q')^B{l,8)-PiQ') (20") 

in the neutron one, to which one must, respectively, either subtract or add the BjSR 
right-hand side: 

C{Q') = 2.IE{Q')-1-C{Q') (21) 

in the first case, and 

cm = 2 • I^iQ') + i • C(Q^) (21') 

in the second. In these cases the last tems are no longer part of the valence dis- 
tribution normalization, and therefore they carry a theoretical error, difficult to 
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quantify , but approximable with that of the most precise experimental determi- 
nation of the BjSR integral (i.e. with that of Ref. 46). 

Taking these aspects into account, the fits yield for the EJSR the integrals 
listed in Table VI for the three valence parametrizations, which show a marked 
reduction with respect the naive expectations for this quantity, tabulated in the 
second column of Table IV. As already said, this is evidence for either presence of 
a sizeable reduction of and/or qq (due to a polarized strange-quark density or 
the breaking of flavour SU(3) symmetry*^^) or an appreciable contribution from the 
axial anomaly (or both). 

From the previous table one can also see that the three parametrizations, 
though behaving quite differently at a; — in the combination g{^^ , do not produce 
neither appreciably different x^'s (the BLP fit turns worse than the other two, wich 
statistically are on the same level, but only by a factor ~ 2.5), nor large variations 
(contrary to naive expectations) for Iq~^"'{Q'^): this is due to the normalization 
of the valence components to the BjSR imposed by Eq. (15)-(16), valid for all 
parametrizations regardless of their behaviour as a; ^ and quite robust, since the 
only isovector contribution from the sea could come, under the hypotheses adopted 
here, from the Pauli principle "penalty factor" . The different behaviours as a; — 
of the valence parametrizations reflect thus only on the quality of the fits, tending 
to prefer distributions which stay finite in this limit, since diverging ones become 
harder to accommodate to the behaviours of gi for neutron and deuteron targets 
at moderate and high values of the Bjorken variable x. 



Table VI 

Integrals of the EJSR from fits to the data 



Expt. 


gVGeV^ 


FRP fit 


SRP fit 


BLP fit 


E142'^ 


2.0 


0.098 ± 0.020 


0.103 ± 0.020 


0.118 ± 0.025 


E143^' 


3.0 


0.073 ± 0.020 


0.068 ± 0.020 


0.068 ± 0.022 


E143'^ 


3.0 


0.094 ± 0.005 


0.092 ± 0.005 


0.096 ± 0.010 


SMC^ 


4.6 


0.103 ± 0.032 


0.095 ± 0.032 


0.086 ± 0.036 


SMCP 


10.0 


0.090 ± 0.032 


0.084 ± 0.028 


0.077 ± 0.028 


SMC^ 


10.0 


0.124 ± 0.024 


0.119 ± 0.023 


0.114 ± 0.029 


EMCP 


10.7 


0.086 ± 0.032 


0.080 ± 0.032 


0.075 ± 0.034 



Any interpretation of these data (beyond the existence of a sizeable, negative 
polarization in the sea) has to make clear the nature of the unitary-singlet piece in 
the EJSR. If one follows the naive parton model picture to the end, and identifies 
go with E, one must also have 7 = from the vanishing of its anomalous dimension, 
so that the EJSR becomes now, since gs = — 3As, 

I^+-{Q') = . Ci(«,) + ^ • Csias)] • S - i • Csias) ■ As + {h.t.)j,^o , (22) 
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which, putting 51 = -{h.t.)i^=o, F{Q^) = | • Ci{as) + ■ Csias) and R{Q^) = 
C8(q!s)/[6 • can be turned into a hnear relation between S and As, 

s = ^'^^^^qI^ + R{Q') ■ A« = MQ') + R{Q') ■ A« , (23) 

to be cvcntuaUy combined with the flavour SU(3) prediction for gs = T, — 3As 
(and with qa = Aw — Ad) to disentangle the quark spin components in the proton 
{modulo flavour SU(3) symmetry violation effects). 

Though simple and attractive, Eq. (22) runs however against the PQCD 
prediction on the evolution of light parton densities^^: when the mass-parameter 
of these latter is neglected, the gluon-to-quark splitting function dictates for qq an 
anomalous evolution and at the same time identifies it with the combination in Eq. 
(5). One is now faced with the additional problem of choosing a normalization scale 
fi^ to relate it to the spin composition of the nucleon. For this, let us consider the 
angular momentum sum rule 

S + 2 • [AG{Q^)+ <L,>] = 1 , (24) 

which has to be satisfied at all momentum scales. Here both AG and < > must 
run with the momentum scale, and from the evolution equation for one can see 
that the first will run like 

AG(g2) = • exp[7(a,(g2)) _ ^{as{^i'))] ■ {AG{fi')+ 

+ ■ S • (exp[7(a,(/.2)) -7(«.(Q^))] - l)} , (25) 

while < L2 >= (1 - E)/2 - AG{Q'^), where the first term (1 - S)/2 is positive in 
all reasonably conceivable models of the nucleon spin structure. 

It is therefore "natural" to assume that < > will vanish at some low mass 
scale, where S < 1 will be balanced by some polarized glue distributed between the 
constituent quarks. If we take as such a scale, then AG{p?') = (1 — S)/2 and 
the EJSR becomes 

- J • Csias) ■ As + ih.t.)i^=o - ■ eMlM - 7(«.(/^'))] • Ci(«.) (26) 

D bit 

The linear relation between E and As, Eq. (23) is still holding, but now F{Q'^) 
is changed including the integrated anomalous dimension factor [1 + 3Q;s(/x^)/47r] • 
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exp[7(Q!s) — 7(q!s(/x^))] to multiply the PQCD coefficient Ci(q!s), and to the correc- 
tion 51 one has to add a positive contribution — equal to minus the last term in Eq. 
(26) — coming again from the anomaly contribution at the normalization scale /U^. 
This last term becomes thus comparable to the EJSR integral Iq~^^{Q'^), so that it 
is not too difficult in this approach to recover a quark content of the nucleon spin 
very close to the flavour SU(3) prediction for ^rg, as people were naively expecting 
before the EMC results^ became public. 

Table VII lists the numerical values for the terms 51 {Q'^), F{Q'^) and R{Q^) 
in Eq. (23) under the two hypotheses: for the second one, /x^ is selected here by 
putting asifJ.'^) = 1, a not unreasonable choice in view of the fact that we expect it 
to lie at the bottom end of the PQCD validity range, due to the smallness of the 
AG involved, whose size at a scale — 1 ~ 2 GeV^ is expected^^ to be comparable 
with < xo >— 1/2. One can note that the PQCD evolution has its largest effect in 
the factor F{Q'^), rescaling Iq~^^{Q'^) to Eq, which in the parton model limit should 
take the value 5/18, while R{Q^) in both cases considered deviates from its parton 
model value 3/5 only at the lowest values of Q^. 



Table VII 

Numerical values for the terms appearing in Eq. (23) 



gVGeV^ 




51 


R{Q^) 




Eq. (22) 


Eq. (26) 


Eq. (22) 


Eq. (26) 


Eq. (22) 


Eq. (26) 


2.0 


0.2310 


0.2315 


0.0145 


0.0509 


0.5566 


0.5553 


3.0 


0.2394 


0.2368 


0.0097 


0.0465 


0.5712 


0.5775 


4.6 


0.2453 


0.2403 


0.0063 


0.0435 


0.5795 


0.5915 


10.0 


0.2520 


0.2442 


0.0029 


0.0404 


0.5874 


0.6062 


10.7 


0.2525 


0.2444 


0.0027 


0.0403 


0.5878 


0.6072 



Using these values one can extract the terms Eo(Q^) in the linear relationship 
of Eq. (23), E = Eo(Q2) + R{Q'^) ■ As, listed in Table VIII. 

To sum the results into single values, the different experiments can be aver- 
aged a la PDG, so as to weight their different (both experimental and theoretical) 
accuracies: the outcomes are 



E = (0.418 ±0.021) + 0.5582 -As (27) 

for the "non-anomalous" treatment of the EJSR, while the correct inclusion of the 
axial anomaly (with the choices ag = 1, < >= at its normalization scale 
Q'^ = ji^) leads to an average 

E = (0.579 ± 0.022) + 0.5796 ■ As . (28) 
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Table VIII 
Values of Sq from the EJSR data 

(averaged over the three fits) 



expt. 


gVGev2 


From Eq. (22) 


From Eq. (26) 


E142" 


2.0 


0.518 ± 0.091 


0.674 ± 0.091 


E143P 


3.0 


0.329 ± 0.060 


0.488 ± 0.061 


£143*^ 


3.0 


0.429 ± 0.025 


0.590 ± 0.026 


SMC^ 


4.6 


0.415 ± 0.134 


0.578 ± 0.137 


SMCP 


10.0 


0.341 ± 0.115 


0.506 ± 0.118 


SMC^ 


10.0 


0.487 ± 0.100 


0.656 ± 0.103 




10.7 


0.331 ± 0.130 


0.495 ± 0.134 



One can note that the treatments proposed here lead in both cases to higher 
values for E than those to be found in the recent literature^'-^^'^^'^°, though of course 
in the second case much higher than in the first. It also noteworthy (in the author's 
opinion) that the second average for Eq is remarkably close to the flavour SU(3) 
value derived above for gs, since estimates on its flavour-symmetry breaking part 
led to expect a value slightly lower than the symmetry prediction^^. 

Note also that this averaging does not produce results much diffierent from 
those obtainable selecting the deuterium data alone, which should be better from 
the point of view of theoretical systematics, since in the handling of the data one 
does not have to correct them with the BjSR integrals. The reduction in the scale 
of the error is of course a product of the use of the full statistics of the (almost) 
complete world data set. 

Intersecting the above error bands on S and As with that from the flavour 
SU(3) result 

£f8 = S - 3 • As = 0.616 ± 0.022 , (29) 
one can now obtain the results {modulo of course SU(3) violations) 

S = 0.373 ± 0.026 and As = -0.081 ± 0.012 (30) 

for the "non-anomalous" treatment of the EJSR, and 

S = 0.570 ± 0.028 and As = -0.015 ± 0.013 , (31) 

for the description of go{Q'^) including the axial anomaly, not very far from the 
naive OZI expectations E ~ gfg and As ~ 0, and practically consistent with them, 
if the theoretical uncertainties^^ entailed by the use of flavour SU(3) symmetry are 
duly considered. 
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Comparing the present analysis with the recently published ones''^^'^^'^^, the 
values of Eq. (31) reproduce within Icr those for the asymptotic limit of the unitary- 
singlet axial charge go{oo)^f'^^ , misnamed E in much of the current literature, 

^0(00)^^=3 ^ _^(^^^^2^)) . p _ 3«s(^^^ _ ^ Q ± 0.025 , (32) 

47r 

but note also that neglecting the anomaly (and trading it for a sizeable As) we 
find a value for the spin sum E still 20% higher than the currently quoted (and, in 
the author's opinion, erroneously derived) values. Note, as well, that, unless one 
assumes SU(3) to be a good symmetry for the axial charges, one can not conclude 
that the low Eq obtained with the identification (70 = E is evidence for a large 
and negative As, since the value 0.418 ± 0.021 lies inside the allowed range of 
flavour-SU(3)-violating effects expected^^ for gg. 

A test of these two possible ways out of the conundrum on the smallness of go, 
in the range of explored by the PDIS experiments, is offered by the possibility of 
measuring the size and sign of AG{Q'^), on which the "anomalous" interpretation 
gives sharp predictions: indeed with the above normalizations one can rewrite Eq. 
(25) as 

AG(Q2) = ■ exp[7(a.(g2)) - 7(«s(/^'))]- 

'^l^^' ^3^) ■ (^^P[^K(/^')) - 7(«.(Q'))] - 1) - ^] } (33) 

which, due to the smallness of the coefl&cient of E for > 2 GeV^, gives quite 
accurate predictions on the gluonic polarization AG in the nucleon, tabulated be- 
low for the value of E given by Eq. (31): these predictions could be tested, for 
instance, measuring the polarization asymmetries for the semi-inclusive process 
I'^N ^ I'^J /tp -\- X, probably a much clearer signature for AG than that proposed 
via semi-inclusive kaon deep-inelastic production^° to test the size and sign of As. 

The third column in this last table displays clearly the effect of the higher 
orders of PQCD in lowering the values of AG at low with respect to the 
leading-order prediction, which would lead to a constant value, even higher than 
the asymptotic limit listed in the last line, due to the large higher-order terms in 
the coefficients Ci^s{oi-s) (see Table II). 
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Table IX 

Values of AG and ctgAG from the anomaly 
contribution to the EJSR and E from Eq. (31) 



Q /Gev 


AG 


ttgAG 


0.215 ± 0.014 


0.215 ± 0.014 


1.5 


0.8453 ± 0.0026 


0.3877 ± 0.0012 


2.0 


0.9996 ± 0.0004 


0.4014 ± 0.0002 


3.0 


1.2025 ± 0.0023 


0.4142 ± 0.0008 


4.6 


1.4030 ± 0.0049 


0.4234 ± 0.0015 


6.0 


1.5224 ± 0.0064 


0.4277 ± 0.0018 


8.0 


1.6843 ± 0.0079 


0.4316 ± 0.0021 


10.0 


1.7439 ± 0.0091 


0.4342 ± 0.0023 


10.7 


1.7725 ± 0.0095 


0.4349 ± 0.0023 


12.0 


1.8208 ± 0.0101 


0.4361 ± 0.0024 


oo {Nf = 3) 


oo 


0.4778 ± 0.0055 



4. Summary and conclusions. 

The several experiments now conducted on the PDIS asymmetries (SLAC 
experiments £142^ and £143^2,13^ ^nd the collaborations EMC^'^ and SMC^^''^^'^^ 
at CERN) do not contradict conventional expectations on the spin structure of 
the nucleon, namely that one should have a strange spin component As much 
smaller than the non-strange-quark sea ones on one side, but on the other no pure 
valence-quark spin components either, as known since more than twenty years (and 
a complete list of references would be almost as long as this paper: it is enough to 
refer to the recent, illuminating papers by Lipkin^^) from the reduction in qa with 
respect to the naive constituent quark model prediction qa — 5/3. 

One finds that it is possible to reduce drastically go{Q^) from the parton 
model plus OZI rule expectation go — gs, due both to the presence of the QCD 
axial anomaly^^ and to its anomalous dimension^^ (and more to the second than to 
the first reason): a correct use of QCD with high orders included (which make these 
effects even larger than the next-to-leading-order alone) is necessary to describe 
the EJSR without confiict with both the experiments and our expectations. The 
"spin crisis" of 1988 was the result as much of an inadequate theoretical description 
as of the low normalization in the F|' values used by the EMC. 

However, the presence of the anomaly is not strictly required by the data^^: 
from a purely statistical point of view, these could be described also by trading it 
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for a sizeable (not necessarily strange) sea component in the nucleon spin: but an 
adequate parton model description excluding the contribution of the anomaly has 
also to exclude the anomalous dimension it carries (for the diagonalization of the 
evolution matrix for the pair S and AG requires one of the eigenvalues to vanish), 
thus increasing by as much as 20% the total spin sum S over currently quoted 
values, and making it possible to explain the data through a sizeable, though not 
impossible, reduction in — = S due to flavour-SU(3)-breaking effects^^. 

Direct measurements of either a large AG or a large As — or even of the 
absence of both, in the case of a large, negative SU(3)-breaking effect — are clearly 
called for to make a definitive choice between these possible ways out of the "nucleon 
spin crisis" . 
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